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^ ' Abstract. The collective-mode spectrum of a multicomponent magnetized ionic 

mixture for small wave number k is studied with the use of magnetohydrodynamics 
and formal kinetic theory. Apart from the usual thermal and diffusive modes, 
the spectrum contains a set of four oscillating modes. By evaluating the 
^3 ■ k^ contributions to the eigenfrequencies, the damping and the dispersion of 

ry>\ ' these oscillating modes are determined. The long-range nature of the Coulomb 

interactions is shown to imply that Burnett terms with higher-order gradients in 
jS I the linear phenomenological laws have to be taken into account in order to obtain 

f-H , a full description of all damping and dispersion effects. 
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1. Introduction 

The dynamical properties of classical Coulomb systems change considerably when 
external magnetic fields are applied. In particular, the spectrum of large-scale 
collective modes is strongly influenced by such fields. Whereas the modes associated to 
the longitudinal and the transverse degrees of freedom may be excited independently 
in an unmagnetized system, these degrees of freedom get coupled when a magnetic 
field is switched on, at least for modes propagating in an oblique direction with respect 
to the magnetic field. As a consequence, both the well-known plasma oscillations and 
the viscous modes get mixed when the isotropy of the system is lost by applying a 
magnetic field. 

The influence of magnetic fields on the mode spectrum of Coulomb systems has 
been studied in several model systems. For the one-component plasma model (OCP), 
which consists of charged particles moving through an inert neutralizing background 
of opposite charge, the changes in the mode spectrum have been studied by means 
of magnetohydrodynamics PP and in the framework of formal kinetic theory PP-PI- 
Since the dynamics of the magnetized OCP model has an anomaly that is caused 
by a divergency in the heat conductivity ^ |S] , a full treatment of the effects of a 
magnetic field on the mode spectrum calls for a somewhat more general model. The 
multicomponent ionic mixture (MIM) in a magnetic field is quite useful in this respect 

u 

In the magnetized MIM model several species of charged particles are present. 
If the charge-mass ratios of the particles are not all equal, the anomalous behaviour 
of the heat transport is no longer present, as has been shown in Ref. ^. In that 
paper the modes of the magnetized MIM have been studied by formal kinetic theory. 
More recently [Jj , a magnetohydrodynamical approach has been used to determine the 
influence of the magnetic field on the modes. The results were found to corroborate 
those found in |B]. In the long- wavelength limit the spectrum of collective modes 
consists of four oscillating modes with a non-vanishing frequency, and of a set of 
coupled heat and diffusion modes with frequencies tending to for vanishing wave 
number. By including terms of higher order in the wave number, the latter modes 
could be discerned. The second-order terms in the mode frequencies describe damping 
that is caused by dissipative transport phenomena. 

In Ref. ini and |7| the mode frequencies of the oscillating modes have been obtained 
in the limit of vanishing wave number only. Hence, dispersion effects for these modes 
could not be studied. Furthermore, damping effects that are of second order in the 
wave number were not treated either. In the present paper we wish to fill this gap 
by evaluating the second-order contributions to the oscillating-mode frequencies in a 
systematic way. 

As in our previous paper jjj we shall use both a magnetohydrodynamical 
approach and formal kinetic theory. Interestingly enough, it turns out that standard 
magnetohydrodynamics is not enough to determine the frequencies of the oscillating 
modes in second order of the wave number. As we will show, so-called Burnett 
terms will be needed to obtain a consistent picture. For systems of neutral particles 
these Burnett terms jSi~ll2j, which describe non-local effects in transport phenomena, 
are of rather limited importance. They play a role only when one is interested in 
contributions to mode frequencies that are of higher than second order in the wave 
number. For charged-particle systems, however, the long-range nature of the Coulomb 
interaction leads to a reshuffling of terms in the magnetohydrodynamical equations. 
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which enhances the relevance of the Burnett contributions. 

The paper is organized as follows. In Section [3 the conservation laws and the 
balance equations of linearized magnetohydrodynamics for the magnetized MIM model 
are formulated. The linear phenomenological laws relating thermodynamic forces 
and flows are presented. In particular, the Burnett contributions will be written 
down explicitly. The Fourier-transformed versions of these magnetohydrodynamical 
equations are given in Section |21 The next two sections contain the analysis that 
leads to our main result, namely the eigenfrequencies of the oscillating modes up 
to second order in the wavelength. In Sections El and \7\ a kinetic approach is used 
to rederive these frequencies and to obtain the kinetic equivalents of the transport 
coefficients, both in leading order and in the higher-order Burnett regime. A few 
auxiliary thermodynamical relations and the basic formalism of the kinetic eigenvalue 
problem are presented in two appendices. 

2. Conservation lavifs, balance equations and phenomenological relations 

An ionic mixture consists of several species of charged particles (with charges of the 
same sign), which move in an inert uniform neutralizing background. If a uniform 
external magnetostatic field is present as well, the motion of each of the particles is 
governed by the electrostatic forces from the other particles and from the background, 
and by the Lorentz force due to the external field. We assume that the motion of the 
particles is so slow that the magnetic fields generated by them can be disregarded. 

In thermal equilibrium the particle density riu of species a (with a — 1, .... s) is 
uniform and time-independent. The total charge density q^ — J^a ^a'^a (with e^ the 
charge of a particle of species cr) is matched by the background charge density —q^. 
The total mass density m^ follows as 'Y^^'m^n^, with m„ the mass of a particle of 
species a. 

If equilibrium is slightly perturbed, the particle densities acquire deviations 
drio-ir, t) that in general are non-uniform and time-dependent. Moreover, each species 
in the perturbed system will get a non- uniform and time-dependent velocity fcr(r, t). 
In the following we assume that the perturbations vary slowly in space and time. The 
characteristic length I of the spatial variations should be large compared to the average 

— 1/3 

interparticle distance Ua- , for all a. The characteristic time for temporal variations 
should be large compared to l/c, with c the velocity of light. 

The continuity equation for species a reads in first order of the perturbations jZj : 

— Sn^{r, t) = -n,V • v{r, t) V • Ja{r, t) (2.1) 

at ma- 

Here v{r,t) = 'Y^^m,anaVa{r^t)/m,y is the barycentric velocity and J^ = 

rUr, na[va{r, t) — v{r, t)] the diffusion flow of species a with respect to the barycentric 

motion. These flows satisfy the identity ^^ Jair, t) = 0. Upon multiplying by e^ and 

summing over a we get the conservation law for the total charge: 

^Sq.ir,t)^-q,V-v{r,t)-J2 (^ ' 7^) ^ ' '^^^''^'^ ^^''^^ 

As an independent set we will choose (j2.2(l in combination with l|2.1|l for cr 7^ 1. 

The time dependence of the barycentric motion is given by the equation of motion, 

which up to first order in the perturbations reads (3: 

d 
my—v{r, t) = -ySp{r, t) - V ■ 6P{r, t) + q^ [E{r, t) + c-^v{r, t) A B] 
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+C-' y (^_^\j^(^r,t)AB (2.3) 

Here dp is the perturbation of the (local) equilibrium pressure, while dP is the 
viscous pressure tensor. The electric field is determined by Sq^ through Gauss's law 
V • E(r,t) = Sqv{r,t). Furthermore, B is the uniform external magnetostatic field. 

The perturbation 6e of the internal energy density satisfies a balance equation, 
which gets the form of a conservation law when only terms up to first order in the 
perturbations are taken into account: 

^fc(r,t) = -KV ■ v{r,t) - V • J,{r,t) (2.4) 

Here h^ = e + p is the enthalpy per volume, with e the (uniform) equilibrium internal 
energy density and p the equilibrium pressure. Furthermore J^ is the heat flow. 

The conservation laws and the equation of motion contain the heat flow Jg, 
the diffusion flows J^ and the viscous pressure tensor 5P. These quantities can be 
expressed in terms of thermodynamic forces by means of phenomenological relations 
|13|. The forces that drive the heat flow and the diffusion flows are [7| 

Xe^y (^) (2.5) 



Xa ^ ^T^ 



V^_ M (£; + !„ A B)(2.6) 
1 \ nicr nil ) c 



Mff Ml 

The subscript t indicates that the gradients of the chemical potentials /io- have to be 
taken at constant temperature T. 

In the standard formalism of non-equilibrium thermodynamics the diffusion flows 
are written as linear combinations of X^ and X„ : 

J, = U, •X,+ Y, Laa'-X,, (2.7) 

'^'(7^1) 

The phenomenological coefficients Lo-e and Lg-o-' are anisotropic tensors of second rank. 
The modified heat flow J'^ = Jg — '^^(h^/m,a)J^, with ha- the partial specific 
enthalpy per particle, is a linear combination of the thermodynamic forces X^ and 
X„: 



J'g - Lee ■X,+ > ; L,a ■ Xa 






(2.8) 


(^(^l) 








As a consequence, the heat fiow J^ itself gets the form 


m 






J, = L„ • Xe + > ; Ua ■ Xa 






(2.9) 


(^(^i) 








with the tensorial coefficients 








L.e - Lee + E f ^^ " ^^ ^ >-., 






(2.10) 


a^l)\^^' ^^J 






(2.11) 


The tensors L satisfy the Onsager relations: 








LlUB)=Lil{~B) , L^(B) = L^(-B) , L^, 


(B)- 


--K'A~B) 


(2.12) 



Oscillating modes of a multicomponent ionic mixture 5 

Similar relations hold for L. 

Finally, in standard non-equilibrium thermodynamics the viscous pressure tensor 
is proportional to the thermodynamic force X^, = —Vv: 

SP = r]:\y (2.13) 

The symbol : denotes a double contraction, so that the ?j-component of the right-hand 
side is 7yy"™X""\ The components of the fourth-rank anisotropic viscosity tensor t] 
satisfy the Onsager relations: 

^ymn(B) = j^mnij,^_Q^ (2.14) 

Furthermore, jf^™-^ is symmetric in ij and in mn. 

For neutral-particle systems the above phenomenological relations are sufficiently 
general to determine the damping and dispersion of the collective modes up to second 
order in the wave number. However, for a system of charged particles this need not 
be true any longer, at least when various species of particles are present. Whereas the 
conventional phenomenological relations are still adequate to determine the modes of 
the OCP model ^"0) the standard formalism has to be generalized for the case of 
the MIM model. For that system the damping and dispersion of the diffusion and 
heat modes up to second order in the wave number are correctly found by inserting 
the above phenomenological relation in the balance equations, as shown in 0. In 
contrast, one has to include Burnett terms with higher-order gradient operators in the 
phenomenological relations in order to determine the damping and dispersion of the 
oscillating modes up to second order in the wave number. 

The relevance of higher-order Burnett terms becomes obvious upon inspecting the 
charge conservation law. When 1)2. 7|l with l|2.6|l is substituted in the last term of 1)2. 2|l . 
the divergence of the diffusion flow J^ leads to a spatial derivative of the electric field. 
The long-range nature of the Coulomb field implies that such a derivative does not 
necessarily yield a contribution of first order in the wavenumber in Fourier language. 
In fact, from Gauss's law it follows that the spatial derivative of the electric field is 
proportional to (5g.y itself, so that it is of the same order in the wave number as the left- 
hand side of (|2.2|) . Hence, if one wishes to determine all contributions to the oscillating- 
mode frequencies up to second order in the wavenumber, one has to go further than 
the standard linear law (|2.7|) . A similar effect is seen in the equation of motion (|2.3|) , 
of which the last term is proportional to the diffusion flows. According to (|2.6|l and 
1)2.7(1 the latter contain electrodynamic terms without gradient operators, so that in 
Fourier language terms appear that are of less than flrst order in the wavenumber. 

From both these arguments it becomes clear that in order to determine the 
damping and dispersion of the collective modes up to second order in the wave number 
one is forced to generalize the phenomenological expression for the diffusion flow by 
including Burnett terms of higher order in the gradients. By invoking parity invariance 
one finds the following general form of the diffusion ffow up to second order in the 
gradients: 

J. = U,.X,+ ^ L...-X.,+ Y. L^^l, ■: x(<;)-HL(<=i i X(^)(2.15) 

a'(#l) <T'(#1) 

with the generalized thermodynamic forces: 

,(c) ^ f ea ei\ ^^,„ , 1 



1 \mcr rrii J c 

Xi^' = - ^VVt; (2.17) 



Oscillating modes of a multicomponent ionic mixture 6 

Both L^'^^, and Lyj are fourth-rank tensors, with the following symmetry properties: 
j^(c^ijmn .^ syjximetric in jm , while L^v™"' is symmetric in its last three indices jmn. 

(c) 

Furthermore, L^^, has the Onsager symmetry: 

^ic)jjmn ^_g^ ^ ^ic)njrm ^_^^ ^2.18) 

The symbol : in (|2.15ll indicates a triple contraction; the i-component of the last term, 

r . , .7- (c)iimn ^t^(c)jmn 

tor mstance, is L^v -^v 

According to H2.15|) with H2.17|l , the diffusion flow depends on the gradient of the 
thermodynamic force X„. Likewise, the viscous pressure gets coupled to the gradient 
of the force X^ when higher-order terms in the gradients are included: 

<5P = ,7 : X„ + ^ Li;) : X^^') (2.19) 

with the generalized thermodynamic force 



j^(c') ^ I f Ba ei \ ^^^ 1 



,ViE + -vAB) (2.20) 

1 yTTicr mi J c 

and a fourth-tank tensorial coefficient LI^J . It is related to LJ,.'^j by the Onsager relation 

Lif^'^'^iB) = L(,<;|"^'™(-B) (2.21) 

which is similar to H2.18|) . As a consequence, Ly^™'"' is symmetric in ijm. 

The generalized phenomenological relations H2.15|l and (|2.19|) contain terms with 
higher-order gradients of the charge density 5qy and the barycentric velocity v. These 
terms represent non-local effects in the phenomenological relations. Clearly, the 
diffusion flow and the viscous pressure at a certain position are partly determined 
by the charge density and the velocity in the neighbourhood of that position. The 
relevance of such non-local effects need not surprise for an ionic mixture in which 
long-range interparticle interactions are present. 

3. Fourier transforms of the balance equations 

To determine the properties of the collective modes we need the Fourier transforms 
of the balance equations. The Fourier transform of the particle conservation law (|2.1|l 
follows upon substitution of H2.15|l . with (|2.5|) . H2.(j|) . H2.16|l and H2.17|l . Up to second 
order in the wave number we get 

7— 5n„ = —ikbjp —^ Jm^ k ■ v — k k ■ L„ ■ [ k — -^ + i 7 ^Jniy v A . 

ot qy mg. T V '^ 



k^ — k-Lar-ksil- 



rUrr \ 1 



-k-^ E fe-U.-fef^-^) (3.1) 

Tin ^ I ' ^ \ Tn _ / 771 ^ l ^ 



nir, 1 ^ — ' V TTlrr' "ll / 'r 

<T'(#1) ^ '" '^T 

Here k = k/k is a unit vector in the direction of the wave vector. Furthermore, we 
introduced the following combination of phenomenological coefficients: 

U = E f — - -) L..' (3.2) 
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The ratio of the plasma frequency ujp — qy/^/rn^ and the cyclotron frequency 
ujc — q^B /{nivc) is denoted as 7 = ujc/ujp- As a consequence of the long range of 
the Coulomb interactions, large-scale charge perturbations are suppressed effectively. 
In fact, one may prove |14) that in thermal equilibrium the Fourier-transformed 
fluctuations Sq^ are proportional to k. For that reason, the 'reduced' charge-density 
perturbation Sqv/k will be treated as being of order fc". It should be noted that upon 
adopting this convention all contributions from the generalized thermodynamic forces 
in (|2.15(l drop out from H3.1(l , as these become of higher than second order in the wave 
number. 

The Fourier transform of the charge conservation law (|2.2|l is derived in a similar 
way. Dividing the conservation law by k so as to obtain the time derivative of 6qy/k, 
we see that the terms with the generalized thermodynamic forces in (|2.15(l contribute 
in second order of the wave number. We get 

d fSqy\ . f ^ f , f f ^1^ , ■ I AD 

+k^-kk : L(^) :fc ( A; ^ + 177^1; A B 

-ifc2^^— fcfc : d^^ fcVWv (3.3) 

y/rriy T 

In the first two lines new combinations of phenomenological coefficients show up: 

L = E Y. (---)(—--) L- (3.4) 

U = E f — - -) L.e (3.5) 

Analogously, the terms in order k^ contain combinations of generalized phenomeno- 
logical coefficients: 

^ .S .^1) v^ ~^i)[^~ ^1) ^^^' ^^-^^ 

Li^) = E (— - -) Li'^i (3.7) 

The Fourier transform of the equation of motion follows upon inserting the Fourier 
versions of (jTT^ and fTT^ in ^^. Using 123-^31, fTT^ - PT7|) and fT^ we 
get 

\/my -—V — —iujy { k —-^ + i 7 ^niy v A B ] +1 — B AL ■ { k — -^ + i 7 Jruy v A B 
at \ k IT \ k 



^-^-{k^BAU-ksll- 



—ikujpk- ikjBALe-kS 



-ifciy BAU-fe ^-^ -k'^k-n-.kv 

T ^ \ ma mi J J, yjmy 
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(c) . f. / f. ^Qv 



k ( k ^^ + ij^my V A B 
k 
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fc-L 

T 

-k^- 



/rriv T 
1 



■BA 



fc-L^,") -.ki^T^v) 



^ k ■ l'^^} -.kik^^ + i7 VWv A B 

m,„ 1 \ k 



with the abbreviation 



E 









(c) 



(3.8) 



(3.9) 



FinaUy, we have to determine the Fourier transform of the energy law (|2.4() . With 
the help of H2.9|l we get 

d 



5e = —ikujp — Jrn^ k ■ v — k — fc-L^- \ k — - + i 7 Jrn^ v A B 
ot qy T \ k 



-rfc-U^ -kS 






'^(7^1) 



(5^CT (5^1 



mi / rp 



(3.10) 



with /i„ = e + p the enthalpy per volume. The tensor L^ is defined in terms of (|2.11|l 
in the same way as in (|3.5|l : 



E 



m„ 



ei 

?7ll 



(3.11) 



The conservation laws (|3.1|l (for cr 7^ 1), (|3.3() . H3.10|l and the balance equation 
H3.8I) form a complete set of coupled linear differential equations. They govern the 
time dependence of the perturbations in the partial densities, the charge density, the 
hydrodynamical velocity and the energy density. In the following, we shall determine 
the oscillating modes of the system for small wave numbers, by solving the associated 
eigenvalue problem in consecutive orders of the wave number. 



4. Oscillating modes in zeroth and first order of the vi^ave number 

In zeroth order of the wave number, the equations l|3.1|) for the perturbations of the 
partial densities n„ and (|3.1U|) for the perturbation of the energy density e get a simple 
form, as all terms at the right-hand side vanish for k tending to 0. Hence, we are left 
with the charge conservation law H3.3|l and the equation of motion (|3.8|l . The latter 
can be decomposed in three equations for independent components of the velocity by 
writing 



v = j^{k\\ ■v)k\i 



fc2 



(fe^ • v) A; J 



1 



[{k±_ AB)-v]kj_AB 



The Fourier variable k has been written as fc 



k±, with B ■ fcj 



(4.1) 



and 



B A fc|| ~ 0. Furthermore, we introduced the notations A;|| = k\\/k^ k± = k±/k, 
fc|| — |fe||| and k± = |fc_L|. 

Using the above decomposition of H3.8() in zeroth order of k and combining the 
resulting three equations with the zeroth-order approximation to (|3.3|l one gets four 
coupled equations from which the oscillating modes 03^ and the associated frequencies 
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■^Ap ' 'with A = ±1 and p — ±1, can be obtained [7|. The time dependence of the modes 
is determined by the eigenvalue equation 

|aS=-i4>S (4.2) 

with frequencies z^J that are solutions of the quartic equation 

>ir - -s - -^^f] [(4°; - i76f - ^'(-^ + ^)'" 

+ 4p'fci {(4°;-i7^'') [5'^-K + 6)2]-2i7K + 6)H'}=0 (4.3) 
The coefficients b, h' and c are defined as 
b=-^L' , 5' = -Il^ , c--- ffc^ill+fclL^) (4.4) 



T ' T ' T V '^" 

with L'I, L^ and L* independent parameters that determine the anisotropic tensor L 
through the decomposition 

L^ L^^ BB + L^ (iJ - Bb) +L^e-B (4.5) 



with e the Levi-Civita tensor. Since b and b' are real, while c is imaginary, the 
four solutions of (|4.3|l come in pairs. Their labels are chosen such that one has 

^Ap — ~^A,-p- 

The four mode amplitudes can be written as [7|: 



"Ap " "IT "^ ^"^" '"^P ' '" 


(4.6) 


The vectors v\p have the form 




'^Ap = 4p ^11 + "Ap ^-L + ^'Ap k±AB 


(4.7) 


where the components are defined as 




II c., ^ 4?(-P + ^) , zi>'-i7[6'2 + (-p + 6)^] 

^Ap (0) ' ^Ap (0) ' ^Ap (0) 
^Ap ^Ap ^Ap 


(4.8) 


with the abbreviation 




AS = (4";-i75f-^'(-^+^)' 


(4.9) 


To determine the first-order corrections to the above oscillating modes 


we write 


«S = «S + fc«S 


(4.10) 



Here and in the following, a superscript '"•' in round brackets denotes the sum of all 
contributions up to and including a given order n in fc. A superscript ["1 in square 
brackets denotes the coefficient of the contribution of a specific order n. In particular, 
one has al = aU . As an Ansatzwe will try and solve a first-order eigenvalue equation 
of the form 

|-a«=-izf 4^) (4.11) 

Qf Ap \p Xp \ J 

with the same zeroth-order frequency as in 14.2|l . and with a first-order term of the 
general form 

«aJ = ^Ap,£ fc + XI ^^".'^ '^"'^ ("^-12) 
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where the coefficients have to be determined as yet. Inserting (|4.1U|) in H4.11|l we find 



d_ (0) 



(0) 



k 



d_ (0) 
dt^'^p 



iW 



d_ [1] 



1(0) 






(4.13) 



up to ffi'st order in the wave number. Of course, the zeroth-order terms drop out in 
view of (|4.2|l . Upon substitution of H4.12|) and inspection of (|3.1|l and (|3.1U|) one infers 
that the last term at the left-hand side vanishes. Hence, we are left with the relation 



dt ^P 



(0) [1] 

Ap <p 



(4.14) 



from which the coefficients Axp j, and ^Ap,o- have to be determined. 

To evaluate the left-hand side of (|4.14l) we need the conservation law for the charge 
density and the balance equation for the three components of the velocity up to first 
order in the wave number. These follow from (|3.3|l and H3.8|l . with H4.1|l inserted. To 
obtain more suitable forms for these equations we rewrite the chemical potentials ^a 
in (|3.3(l and H3.8|l in terms of new chemical potentials jia = fJ-a — (ecr/ei)/ii jl5| . These 
are better suited to the description of an ionic mixture, which satisfies the constraint 
of overall charge neutrality. The resulting expressions for the perturbations of the 
chemical potentials are given in ljA.2|) . Furthermore, we decompose the tensors L^ 
and Lg. in a similar way as in (|4.5(l . In this way we arrive at the following forms of 
the balance equations: 



d_ 
dt 



Sqv 
k 



.d_ 
dt 



k 



fell • t) — i (ujp + b) y/m.v fej 



-i b' y/m^ {k± A B) ■ V 



-ifc 



Tc,S 



T 



T J2 c„6 

<^(#i) 



(fe||.^) 



> . f2 Sqv 

V I = -lWpfe|| — 



fc^^ 



^pft,ll 



qv 



Sp 
qv 



(4.15) 
(4.16) 



d 



Sqv 



fmv -^ f fe± ■v\ = -\ (cjp + b) fci ^ + 7 6' 



'k± ■ V 



—7 {ujp + b) ^Jmv (fe_L /\ B) ■ V 



+ifcfc^ 



Tb^S 



'^(7^1) 



{ujp- 



5)^ 
qv 



(4.17) 



" dt 



(fcj 



AB 



ib'k 



! £.2 ^ 
k 



+ 7 {ujp + b) y/m^k±^ 



-7 b' ^JmZ {k± AB) -v 



-i k k", 



Tb'S 



T 



T Y.KS 



b'^-p 

qv 



(4.18) 



The abbreviations 6^, 6^, b^, b'^, c^ and c„ are combinations of the components of L, 
Lp and L^- 7 : 



T 



■Ll 



61 



qv 



T 



b' ■ 



(4.19) 
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6' = 



7 



H 



kWl + fc^ Li 



Qv 



niaT 



■m„T 



L 






kILi 



Qv 



11 
(4.20) 

(4.21) 



with (7^1. 

The terms of order k in (|4.15() - (|4.18|l contain the perturbations S{l/T), 5[jla/T) 
and 5p. These can be expressed in terms of the perturbations 5n„ and 5e by employing 
the thermodynamic relations l|A.9|) . IJA.10|I and ljA.13|) of [Appendix A| of which only 
the terms of order fc° are needed for the moment. Substituting the results in (|4.14|l 
with (|4.t)|) . we arrive at an expression for a}-^'^ which is indeed of the form of (|4.12|) . 
Upon using (|4.3|1 we obtain the following expression for the coefficients: 



A 



^p,j 



r' 



Wqv 



m^-s, + J2m-,}f,,, 



(4.22) 



with j and j' equal to e or a-{^ 1), with /3 = l/{kBT), and with M--, and Sj defined 
in (|A.4|1 and ljA.8|l . Furthermore, we introduced the abbreviation 

1 



-^Apj - 



,(0) 



Ci + Hpbj-vipb'^)kl 



with components v\ defined in 



(4.23) 



5. Oscillating-mode frequencies in second order of the v^rave number 



Having determined the oscillating modes up to first order in the wave number, we 
now turn to the contributions of second order in k. In particular, we want to evaluate 
the mode frequencies in second order, as these determine the dispersion and damping 
effects. The mode amplitudes in second order are less interesting. Fortunately, their 
explicit form is not needed in the following, as we shall see. 

[21 

The second-order contributions z^' to the eigenfrequencies of the oscillating 
modes follow from the eigenvalue equation 



d_ (2) _ _ 

at "^p ~ ' 

with the mode frequencies 

J2) _ (0) 
^Ap ~ ^Ap 

and the mode amplitudes 



• (2) (2) 
-l^p ^Ap 



'Xp 



(5.1) 



(5.2) 



^Ap - 



^Ap 



^Ap 



,[2] 



(5-3) 

The zeroth- and first-order terms in the amplitudes have been given in (|4.6|) and (|4.12|l 
with (|4.22|) . The second-order terms have the general form 

a'-^l = Axp.q—^ + A\p^/mZk\\ ■ v + Aj^^ ^/rri:^ k ±_ ■ v + A\p^/mZ{k±_ A B) ■ v 

with coefficients A^p that are as yet unknown. 

Substitution of H5.2|) and H5.3|) into (|5.1|) yields the second-order equation 



(5.4) 



[ ^ „(o)l 


[2] 


P [1]] 


[1] 


\ 9 [2]1 


[si"^pj 


+ 


[dt "H 


-1- 


i^^H 



(0) 



• (0) [2] 
-l^p <p 



[2] (0) 
4p «Ap 



(5.5) 
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The right-hand side is a Hnear combination of the four independent perturbations 



Sqv/k, 



k 



/rriy 



k± ■ V and ^TOi, {k± AB)-v, as follows by inserting H4.6|l - H4.7|l 



and H5.4|l . Likewise, the third term at the left-hand side of H5.5(l can be written in terms 
of these independent perturbations by inserting H5.4|) and substituting the zeroth-order 
contributions of H4.15() - (|4.18l) . In the second term at the left-hand side of (|5.5|l we 
substitute (|4.12|l and use (|3.1|l . (|3.10|l . The first-order contributions in the latter 
equations can be rewritten as linear combinations of the independent perturbations 
in the way shown in (6.1)-(6.2) of [3- Finally, we have to consider the first term at 
the left-hand side of H5.5|l . Formally, it can be written as a linear combination of the 
independent perturbations as well: 



d_ (0)' 



D 



\p,q ■ 



SQv 



D 



Ap' 



k\\-v + D 



Ap 



fci-t;+L»L\M^(fc_LAB)-t;(5.6) 



The coefficients D\p will be determined later. 

Having taken the above steps we are in a position to compare the contributions 
of the four independent perturbations at both sides of H5.5(l . Upon equating the 
coefficients we arrive at the following set of equalities: 



(5.7) 



ic-zf^)A,p.,+copklA 


yip + [ujp + b) fci 


Aip-b'klAip 


= 4'j-i^Ap,, 


-Aap,£ Ce + C - > ^Ap.cr C„ + C 


A (0) .11 [21 II . r^ll A ^-" Sr^ A '^o- 

UJp Axp,q - z^ A^ = z^ v^ - 1 £'ap - ^Ap.s '^P~- y , ^^p-" ^P — 


(cp + h) A^p^, + iijb'~ z^^} )Aip + \^ {ojp + b) A\p = 7^11 vip i Dip 


-Ap,e 


b, + {iUp+b)^ 




6, + (wp -f 6) ^ 


6' Aap,, - i 7 (^p + 6) Aip + {nb'~ zi°; ) A{p ^ zgj v{p - i Dip 


-Axp,e 




r Axp.a I b'^ + 





(5.8) 



(5.9) 



(5.10) 



The above system of equations for the four unknowns A\pq^ ^a«' ^Ap ^^"^ ^\p 
turns out to be a dependent set. In fact, the determinant of the matrix of coefficients 
appearing at the left-hand sides is exactly equal to the left-hand side of the quartic 
equation H4.3|) . so that it vanishes. The mutual dependence of the equations is a 
mathematical reflection of the fact that the amplitudes of the modes are fixed up to 
a multiplicative constant only. Because the system is a dependent set, the matrix has 
got a left-eigenvector with eigenvalue 0. This four-dimensional vector is found as 



1 



Ap II 



"Ap '^i J 
J2] 



-^Ap '^± 



(5.11) 



The second-order contributions zP to the eigenfrequencies of the oscillating modes 
now follow by taking the inner product of (|5.11l) with the set of equation (|5.7|l - l|5.10|) . 
The left-hand side of the resulting equation vanishes by construction, so that we are 
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left with an equality from which zP can be determined. Using l|4.3|l to eliminate c we 
find: 



'Ap 



Nx 



P \ ^Ap 



qvj ^^_^^^ V qv 



+iDxp,q + ^v\pklDl^ + xv 



\p ^± D\p i ^\p ^± D\p 



(5.12) 



where the constant N\p is defined as 



N: 



\p 



+ K)'^^Kp)'fci-Kp)'fci 



(5.13) 



The expression (|5.12() still contains the coefficients Dxp that have been introduced 
in 1)5. 6|l . These will now be determined by substituting (|4.6|l in the left-hand side of 
(|5.6|l . Employing the balance equations (|4.15(l ~ (|4.18|l and using the thermodynamical 
relations of [Appendix A| we encounter contributions that follow from the terms 
proportional to 6a., in (|A.9|l . (jA.10|) and (|A.13|) . The ensuing term in i?Ap,9 is: 



^^p,q ~ Q ^Xp 



E^v,.^. 



2/3g. 



■W, 



2(7„ei 



(5.14) 



A further contribution D'^ to D\p^q follows by considering the terms of order k^ in 
(ins and lESl). These yield: 



n" = _ 



' vxpk : \S;;} ■.kk + ^(^k-i-fVxpAB]k: L^^^ : kk 



Likewise, we find for the other coefficients: 
1 ." .".'■ i 



D 



\p 



K = - 



my k 
1 



2 vxpk : T] : kk\\ - 



vxpk : Tj : kk± 



/ TTii) J. /Cii 
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(k ~ i-fvxp ABJk: L^''^ : A;fc|| 



(5.15) 



(5.16) 



rfc2 



k-i"/vxpAB] k : L^"' : k{kj_ A B) 



(c) 



/rn^Tk^ 

7 
/m,y Tk\ 



2 ^k~i-ivxpABjk: L^'^ -.kki^ 
vxpk : l|,^^ : k{k^ A B) 



(5.17) 



^\p 



TTiy fC I 



^ Vxpk : T] : k{kj_ A B) 



17 
Tkl 



k-i'-fVxpAB)k: L'^"^ : kk± 



/rriy Tk'^ 

7 



(k - i-/vxp Ab) k: L^'^) : k{kj_ A B) 
vxpk : L^*;^ : fcfc_L 



(5.18) 



^rriv Tk'j^ 

The fourth-rank tensors in these expressions may be decomposed in terms of their 
invariant parts in a way analogous to H4.5|l . The expressions for the coefficients Dxp 
that result by substituting these decompositions are rather unwieldy and will not be 
presented here. 
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Apart from the coefficients Dxp, the expression H5.12|l also depends on the 
coefficients A\p for which the exphcit form (|4.22|) is available. Upon inserting the 
latter and combining the resulting contributions with those from H5.14|l we arrive at 
the final form for the second-order terms in the eigenfrequencies of the oscillating 
modes: 






C's. 2^f 



„2 



E(:7|r^^e7--^.Vv. + ^E^v,, 



Loy P^ \2Pq, '= V ''' P 



M7,, F: 



^p,j' 



Ap.q + i Ap ^\ Dip + i vip kl Dip - i v\p ki ij^p 



-\D"+W{k\D\^ + \vipk\Dip-\v\pk\D\p\ (5.19) 



where we have used ljA.18|l as well. Our result shows that the damping and the 
dispersion of the oscillating modes in second order of the wave number is caused 
by several rather different mechanisms. The first term, with the sound velocity Cg, 
depends on equilibrium properties only. Of course, dispersion of collective modes 
that is governed by the sound velocity is a well-known effect, both in plasmas and 
in systems of neutral particles. The remaining terms in the first line are linear or 
quadratic in F\p. According to (|4.23|) these coefficients are linear combinations of 6j, 
6' and Cj, which have been defined in H4.19|l -- (|4.21|) with 14.4|l . They are governed by 
non-equilibrium processes through the components of the phenomenological tensors 
L, Lg and L(j, which are linear combinations of \-aai and \-ea- The latter occur in 
the linear laws H2.7|) and H2.9|) for the diffusion flows and the heat flow. Finally, the 
second line of (|5.19|) depends on non-equilibrium processes as well via the coefficients 
D\p. From H5.15|l - I|5.18|) we infer that viscous effects play a role here. However, 
these are not the only dissipative phenomena that are relevant for the damping and 
the dispersion. Higher-order couplings between thermodynamic forces and flows, as 
described by Burnett terms, are important as well. 

For the one-component plasma the expression H5.19|l gets a simpler form. For that 
system -F\p.j vanishes, so that in the first line only the term with the sound velocity 
survives. Furthermore, the Burnett contributions disappear from the coefficients D\pi 
so that only the viscous terms in (|5.16() - (|5.18() remain. The resulting expression agrees 
with that found in pP and |2]. 

6. Relation to kinetic theory 

The collective modes of an ionic mixture can also be derived with the help of formal 
kinetic theory jSJ [7j . In that approach the starting point is the set of microscopic 
balance equations. From these one obtains the collective modes by employing a 
projection-operator technique. 

The microscopic partial particle density is n^{r) — ^^^{r — r^a), with Taa 
the position of particle a of component a. The Fourier transform (5n™(fc) of the 
perturbation n™(r) — (n™(r)) of the microscopic partial particle density , with (n™(r)) 
its uniform thermal equilibrium average, reads 

5n-(fe) = ^e-''='--° (6.1) 

The Fourier transform of the perturbation of the microscopic charge density is 
b<f^(}i) = '^a- e'ijSn';^{k). Furthermore, the microscopic partial momentum density 
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of species a in Fourier space is given by 

g-(fc) = 5]p.„e-'^'-'° (6.2) 

a 

with paa the momentum of particle a of component a. FinaUy, the perturbation of 
the microscopic energy density in Fourier space wiU be denoted by (5e™(fc). 

The conservation laws for Sn^{k) and (5e™(fc), and the balance equation for the 
microscopic total momentum density g"^{k) — J2a9cr^i^) read [HI: 

CSn'^{k) = -—-g^{k) (6.3) 

rricr 

Cg^^ik) = k ■ dr'^ik) - q.,^5C{k) - i V ^^g^(fc) A B (6.4) 

(7 

CSe'^{k)^-k-j^{k) (6.5) 

with (5T™(fe) the Fourier-transformed perturbation of the microscopic pressure tensor 
and J™(fc) the microscopic energy-current density in Fourier space. The Liouville 
operator in phase space £ follows by writing the time derivative F of an arbitrary 
function F as -F = iCF. 

The collective modes are specific independent linear combinations of the 
microscopic basic quantities 

a^(fc)e |^<Sgjr(fe),g"^(fe),&-(fc),<5n-(fc),...,,5n-(fe)| (6.6) 

with i — 0, 1, 2, 3, e, cr (t^ 1). The microscopic charge density is divided by the wave 
number, since the fluctuations in the charge density vanish in the long-wavelength 
limit, as has been remarked already in Sectional Indeed, in leading order of the wave 
number one has 

i([<5g,-(fc)]*<5g-(fc)) = ^ (6.7) 

with V the volume of the system. The fluctuation formulas for the ionic mixture have 
been given in |H] and ^^. 

To derive the collective modes one may use a projection-operator technique. The 
one-sided Fourier transform of the time-dependent collective mode a™(A;,t) is given 

by 

aT{k,z)^~i dte-'aT{k,t) = ^aT{k) (6.8) 

Jo Z + L, 

which is regular for z in the upper halfplane. By introducing a projection operator P 
that projects an arbitrary phase function f{k) on the space spanned by the basis H6.6|l 
and its complement Q = I — P, one arrives at an eigenvalue problem for the modes 
a™(fc) and their frequencies z™(fe): 

^zf^ik)^{[a^{k)Y aTik)) (6.9) 

for all basis functions a^{k). Near z = 2™(fc) the projection of 1)6. 8|l has a pole 
structure of the form 

P -^ aT(k) ^ \—- aTik) (6.10) 

z + L z — z^[k) 
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The eigenvalue problem Ht).9|) leads to a set of s + 4 modes. Apart from a heat 
mode and s — 1 diffusion modes one finds four modes that are the kinetic counterparts 
of the oscillating modes. In zeroth order of the wave number their frequencies z^^ 
follow from a relation of the form H4.3|l . The parameters 6, b' and c are replaced by 
z-dependent functions which follow from the identities \7\ 



k 



CO i Q^ ^C(fc) .(o) ^ .. 



(3 



/rriy V 



{[g^{k)r£Q- 



1 



{[g"\k)rCQ- 



1 



Q£Q^^^^>*"^ - b{z) k± + h'{z) k^AB 
\(o) =_ 



(6.11) 
(6.12) 



QCg''\k)Y^> = -i-fb[z) e B + ijb'{z) (U - BB) (6.13) 



m^F"'' '""" '^ z + QCQ 

The consistency of the second and third of these definitions has been demonstrated in 



The amplitudes a™ (fc) of the associated modes in zeroth order of the wave number 

are given by 

Sq^ik) , 1 



m(0) 



(fe) 



vZ{k)-g^{k) 



(6.14) 



The vector d™ (A;) has the same form as (|4.7|l ~ (|4.8|l . with the parameters b, b' and c 



,m(0) 
-Ap 



replaced by the corresponding functions b{z), b'{z) and c{z) for z ~ 

The first-order contributions to the mode amplitudes have the general form (cf. 



'"W^^^-A- .fc-(fc)+ J2 A'j:,,Jn'^{k) 



"Ap 



\k) 



Xp,e 



(6.15) 



•^(^i) 



with coefficients A™ ■ that have to be determined by substituting H6.14|l and (|6.15|l 
into (|6.9|) for i = Xp and j = e, cr. Employing the fiuctuation formulas |S], ^ and 
the thermodynamic formulas of [Appendix A| we find a set of inhomogeneous linear 



equations for the unknown A™ , with a coefficient matrix given by ljA.4|) . Solving 
these equations we arrive at expressions for A™ ■ of the same form as (|4.22|l , with 
Fxpj replaced by F^nj, with a similar structure as in H4.23|l . The kinetic counterparts 
of (|tTOI) - (|i:^ are 0: 



V 



<5C(fc) 






k 
sc{k) 



CQ 
CQ 



1 



P 



Irtiv V 



{[g"\k)rCQ 



z + QCQ 

1 
z + QCQ 
1 



Q£ 
QC 



fc'"(fc) - —Sq'^ik) 
Sn'^ik) - -^6q^{k) 



)W=c.(^) 



:QC 



([g-(fc)]*/:g 



z + QCQ 

b^{z) kj_ + b'^{z) k± /\B 
1 



fem(fc) _ ^Jg-(fc) 



qv 



:QC 



z + QCQ 

= b,{z)k^ + b'„{z)k^hB 
1(1) 



5n-;{k) —5q^{k) 



(6.16) 
(6.17) 

(6.18) 
(6.19) 



The oscillating mode amplitudes a^^ up to first order in the wave number have 
been obtained as the sum of 16.14|l and l|6.15|l . In the following we shall also need the 
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thermal and diffusive mode amplitudes up to first order in k. These have been derived 
in [7] as 






1(1) 



(k) ^ Sn'^ik) - — SCik) + k 



1 






qv 






The coefficients A™ have the decomposition 



with 



A" 



okf, \1 



Af = 


A-l 


fc| 


+ Af^ k^ - 


b'/kl 


~Cj 


) 


, Af^^ 


^f = 


ib 

7 


j{uj 


p + b) + b'^b' 


{Ujp 


+ 6)2 + 6'2 



A^-^* k^AB 



i bjb' -b'j{ujp + b) 



(6.20) 
(6.21) 

(6.22) 



(6.23) 



for j = £, cr(7^ 1). The functions b, b' , bj, b'^ and Cj are given by (|^Tl|l . l|^TTC|l - l|0^ . 

Here z is equal to the infinitesimal value iO, since the frequencies z" of the heat and 
diffusion modes vanish in lowest order of k. 



7. Kinetic derivation of the frequencies in second order of the wave 
number 

The frequencies z™ in second order of the wave number can be obtained without 
having to solve the full eigenvalue problem (|6.9|l in second order. In fact, it can be 
shown (see [Appendix B| | that knowledge of the mode amplitudes up to first order is 
enough to evaluate the second-order frequencies, since the latter follow as: 



i[2] 



1(0) 



V 



1(1) 



(fe) 



<'''\kr'~^i 



1(1) 



(fc) 



Caf'\k))^'^ 



1 



1(1), 



(fc) CQ- 



1 



m(0) 



QCQ 



-QCaf''>{k)) 



Wi 



(7.1) 



In this expression we have introduced the adjoint modes a™(fe). These adjoints are 
specific independent linear combinations of the basis functions a°{k) 1)6. 6|l . They are 
defined by the orthonormality relations 

^{[aT{k)raf{k))=5,, (7.2) 

From equation \I.\\i it follows that the adjoints of the modes H6.14|l - H6.15|l up to 
first order in the wave number have to be determined. In deriving these adjoints from 
H7.2|) it will be assumed that the transport properties b{z), b'(z), bj{z), b'Az) and Cj{z) 
are slowly varying functions of z, which can be treated as constants. In zeroth order 
the adjoints a™ (fe) are given by the general expression [S] 



^(0)^1.^ - 



a-;'^(fe)^/3(iV-)' 



SCik) 



/mv 



vUk)-g^{k) 



where the vector u™ (fc) can be written as 



^Xp 



^xp '^W + ^\p 



k± + vTp k±/\B 



(7.3) 



(7.4) 
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The normalization constant N^^ and the coefRcients w^J , v'^J- and v™^ follow from 
the orthonormality relations in zeroth order: 

^([arr(fcr<fHfc))(")=<5.A. (7.5) 

with X — ±1, p — ±1, jj, — ±1 and a = ±1. With the help of the fluctuation formulas 
in] these orthonormality relations can be brought into the form 

^a";, {1 + K;(fcr • <.(fc)} = s^,s,, (7.6) 

The coefRcients u^ , v^J- and u™* may be obtained from this equality by choosing 
combinations of indices (Ap) that differ from (p-cr). In that case the orthonormality 
relations read 

(^^5'r<i'^| + «p^)*<fci + «r<^i = -i > (Ap)^(A^c7) (7.7) 

Let us now consider two different solutions z™ and z™^ of the quartic equation 
H4.3(l . An alternative form of the latter is 

c - z';:^ = -c.,<jlfc^ - (.;, + foXfci + b'vTfkl (7.8) 

Subtraction of this equality from its analogue with (Ap) — > (p-a) and division by 
z'^p ' — z"a ' yields the identity 

«r>"i'fc' + <<fci-<<fci = -l , (Ap)7^(M^) (7.9) 
By comparison with (|7.7|l one infers that the coefRcients v'^^ , v™^ and u™* are given 

by 

vTi^KpT , < = «r , < = -«)* (7.10) 

The normalization constant N™ follows from 17.6|l for (Ap) ~ {lJ-(y)- It has the same 
structure as (|5.13|1 . 

The contribution of first order in the adjoint a™ (k) takes the general form 



a^j,'\k)^P{NT,) 



i^,^, fc-(fe) + J2 K.'^ ^<ik) (7.11) 

The coefRcients Af ^ and A™ ^ can be obtained from the orthogonality relations 

^([«rr(fc)]*a7<^Hfc))«=0 (7.12) 

with A == ±1, p = ±1 and j = £,(j{^ 1). Substituting ^^, ((7111) and ^^, 
H6.21|l into H7.12() and employing the fluctuation formulas ^, J4, we get a set 
of inhomogeneous linear equations for the unknown coefRcients A™ . Using the 

expressions (|6.23|l for the coefRcients A™ , A"^^ and A'°* we obtain the solutions: 

K. - (^ApJ* (7-13) 

with j = e,a{^ 1). 

Having succeeded in deriving the adjoints up to first order in k we are ready 
to evaluate the terms of second order z™ in the mode frequencies. In the first 
term of (|7.1|) we employ the fluctuation formulas and the thermodynamic relations of 
[Appendix A| up to second order in k. In particular, we use the Stillinger-Lovett relation 
for the ionic mixture, as given in formula (6.6) of (14) . and the closely connected 
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fluctuation formula involving the pressure and the charge density (sec (7.9) of |14,). 
As a result we get from the first term of H7.1|l : 



N 



in(0) 



4£ 



\p ^\p 



^f3 Y^ A',-^,^. M,,, A'^^,^., - 2 Y^ A'^^,^. M,,, s, - 2A^^^, — 



3,3 



JJ' 



3(j„ 



.(,^1) ^" ^ 3:3' ^" ^^^ 



(7.14) 



A more useful form, which partly resembles (|5.12|) . is found by employing (jA.ll|l and 
the kinetic analogue of 1)4.22(1 : 



N^l z 



m(0) 



Xp '^Xp 



qvj ^^_^^^ \ qv 



+UT.P^.W3 



/3 



2/391, 



W. 



2qvei 



(7.15) 



In the second line we recognize the kinetic counterpart D™ ' of H5.14|l . 

In the second term of 1(7. l|l we insert the laws ((f).3|l and ((6.4(1 . Using the fluctuation 
formulas once again, we find: 



^N^p^P 



■A 



Xp.,£ 



E?^>.. fe- 



"^(^i; 



Qv 



Xp 



(7.16) 



Finally, we have to evaluate the last term of ((7.1(1 . Inserting ((6.14(1 . 1(6.15(1 . 1(7.3(1 
and ((7.111) . we obtain a sum of several contributions. The first of these is obtained 
upon substituting the zeroth-order mode amplitudes ((6.14(1 and ((7.3() in 1(7.1(1 . Using 
(|n31)-(E3I) and the identities QSq^{k) = and Qg"^{k) = 0, we get 



^iV,- (fc-i7t;^;AB).( 



m(c) 



=)■(' 



1 



T^/rrTy 
1 



fc-L:"(^).fc).i;^7 



i 



k-V^^"' -k) ■ [k + i-fv'J^; ^B 
NT,{k~nvl^^AB)-[ 

NTp <p • (k ■ C'^) . fe) . (fe + i 7 vZ* A B 

k-v'^-kj-v^; 



^^Xp^Xp 



(7.17) 



Here we introduced the microscopic equivalents of the transport quantities \-'^l,^ L^'^J, 
L„^' ((7, tr' 7^ 1) and -q. They are defined as 



ym-im Q.^Q^Q 



' Q{g^{k)y)^'^=ikJk-C:^';\z)-kY 



1 

V 

1 

V 

I 

V 



{ {k ■ (5r'"(fc))^' 

([(fc-(5T'"(fc))^' 



^ Q(fe.r-(fc))^)W=fcBffc-CJ^)(z)-fc)'' 



z + QCQ 



Q 



1 
z + QCQ 



Q(fc-x'"(fc))^)(°) = -i(fc-r;"^(z)-fcy 



(7.18) 
(7.19) 
(7.20) 
(7.21) 
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with i, j = 1,2,3. In (|7.17|) we employed reduced transport quantities that are defined 
in a way analogous to H3.6|l . H3.7|l and H3.9|l . All these quantities are to be taken 
The contractions of the vectors k in H7.17|l - H7.21|l arc to be taken 



a.t z — z 



m(0) 
Xp 



with respect to the two inner indices of the fourth-rank tensors; the right-hand side 
of H7.18|l . for instance, equals ifcsfcP (L™j:, (z))*^''-' fc''. The expression H7.17|l can be 
rewritten by inserting 17.4|l with (|7.1()|) . In this way we get: 

„.mt £.2 



^N'^piD 



Xp,q 



where the coefficients D^^„ " 
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^Xp '^± ^Xp 
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(7.22) 



'Xp 

-D"„ , D™J^ and Z?™^ have the same form as their 



magnetohydrodynamical counterparts (|5.15|) - (|5.18|) . 

Further contributions from the last term of (|7.1|l arise by inserting either (|5.14|l 
and fTTT^ or (jOB|) and ^^. Employing (jOT|) - (|Fl^ and (JH1B|| - |H1^ we get: 
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(7.23) 
where H7.1()ll and H7.i;-{|l have been used. As before, all quantities b, b', c, bj, b'j and Cj 



are taken for z — 
in terms of F- 
As a result 



z^^ ' . By means of the kinetic analogue of (|4.23l) we may write Cj 
Furthermore, the quartic equation H4.3|) can be used to eliminate c. 
gets the form: 

If"}} \ flu 
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(7.24) 



By adding H7.15|l . H7.16|l . H7.22|l and (|7.24|) we find an expression for the second- 
order frequencies that is very similar to (|5.12|l 
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(7.25) 



Note that in this kinetic derivation the coefficients A^„ ■ and F?\ ■ take the role of their 

'^P^J ^P^3 

magnetohydrodynamical counterparts ^Apj and -Fapj (with j = e, a{^ 1)). Likewise, 
the coefficients D'^ , -D^ , D™^ and D™* replace their magnetohydrodynamical 
equivalents H5.14|l - lj5.18|l . Upon substituting the kinetic analogues of 14.22|l in (|7.25|l 
we finally arrive at the kinetic expression for the second-order contributions to the 
eigenfrequencies of the oscillating modes: 
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(7.26) 
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This expression has the same form as (|5.19|l . Hence, the kinetic approach corroborates 
the magnetohydrodynamical resuhs for the damping and the dispersion of the 
oscillating modes, at least up to second order in the wavelength. 

8. Conclusion 

In this paper we have shown how the spectrum of collective modes of an ionic 
mixture with an arbitrary number of species can be analyzed by means of 
magnetohydrodynamical and kinetic methods. We have derived explicit expressions 
for the eigenfrequencies of the oscillating modes up to second order. Since the second- 
order contributions to the frequencies determine the damping and the dispersion of 
the collective modes, our results give useful information on the dynamical behaviour 
of ionic mixtures. As we have seen, the damping and dispersion phenomena are 
governed by dissipative mechanisms that go beyond those of standard non-equilibrium 
thermodynamics. In fact, the well-known phenomenological laws for the diffusive flows 
and the viscous pressure are not sufficient to account for the full effects of dispersion 
and damping of the oscillating modes. Whereas the dissipation implied by these 
standard laws yields the damping and dispersion of the heat mode and the diffusion 
modes, as shown in our previous paper 0, we had to generalize the standard laws 
by including higher-order Burnett terms in order to get a complete picture of the 
damping and the dispersion of the oscillating modes up to second-order in the wave 
number. 

Burnett terms have been discussed previously in the context of systems of neutral 
particles |H1-C21- For such systems the Burnett terms generally lead to small effects 
in the dynamical behaviour, so that they can usually be neglected. In the present 
case of an ionic mixture that is no longer true: in a systematic treatment of the 
collective modes the Burnett terms can not be disregarded. The technical reason for 
the necessity to incorporate Burnett terms when investigating the dynamics of an 
ionic mixture turns out to be the long-range nature of the Coulomb interactions. In 
an ionic mixture the long-range electric fields generated by charge fluctuations are 
driving agents in the thcrmodynamical forces. These electric fields lead to terms of 
zeroth order in the wavenumber after Fourier transformation. To be consistent up to 
second order in the wavenumber one has to include terms with gradients of the electric 
field in the thermodynamic forces. Hence the phenomenological laws, which connect 
these thermodynamical forces to the thermodynamic flows, must be used in a form 
which accounts for such field gradient terms. We have to conclude that Burnett terms 
can not be neglected in systems with long-range forces. 
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An ionic mixture satisfies the global constraint of charge neutrality. Therefore, it is 
convenient to introduce instead of Ha a set of new chemical potentials: 

/iCT=M<T -f^i (A.l) 

ei 

for cr ^ 1. As shown in |15| . these potentials arise naturally as Lagrange multipliers 
in a grand-canonical ensemble with a charge- neutrality constraint. The combination 
of perturbations of the chemical potentials fi^ occurring in the conservation laws and 
the balance equations of section |31 can be rewritten in terms of these new chemical 
potentials as [7]: 
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(A.2) 



for a^l. 

The perturbations S{l/T) and 5{jl„, /T) (with a' ^ 1) in (|A.2|I have to be 
represented as linear combinations of the basic perturbations 5e, Sn^ and Sqv This 
can be achieved by starting from the identity 
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with the matrix: 

-de/df3 -dug./dp 
-dfig/dp drig / d{l3ila') 

and the vector 

de/dq^ 
dria/dqv 

All partial derivatives are defined in terms of the independent variables /3 — l/{kBT), 
Pjlrr (with cr 7^ 1) and q^. We used the Maxwell relation de/d{(3fla) = —dn^/df]. 

The elements of V can be evaluated by substituting the equation of state 
p = ^e + ^nkgT into the formulas (3.22) and (3.18) of rcf. jT3], and replacing n 
by —J2a(^i)(^<^/^^ — 1) ?^cr + qv/^i- Treating /?, /3/icr (c ^ 1) and q^ as independent 
variables, we get 
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(A.7) 

(A.8) 
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Inverting HA.3|) and substituting HA.6() and (|A.7|I we finally get 

k 



6(3 = M-^He + ^ M;-„^5n^ - kwj-^ (A.9) 



K#i) 



(5(/3/i,) = M-/fc + Y. Maa'^^'y-kwJ-^ (A. 10) 



^'(#1) 



with Wj defined as 



r ^" <T'(#i) ^^ 

for J and j' equal to e or ct(7^ 1). 

The variation in the pressure 5p can likewise be expressed in terms of variations 
Se, 6n„ and 5q^/k. In fact, the perturbation of the pressure follows from the equation 
of state as: 

Upon inserting the expression (|A.9|I for 5(5 we obtain 

1 _ - H^^% (A.13) 

The inverse matrix element M^^ is directly related to the specific heat Cy per 
particle: 

c. = rf^) =^=lkB (A.14) 

with s the entropy per particle. Likewise, one may derive that the specific heat Cp per 
particle at constant pressure is given as 

Furthermore, the isothermal compressibility kt can be introduced by any of the 
following equivalent definitions: 

1 (dqy\ 1 (dn\ 1 (dmy\ 

where use has been made of the identity dn^ — (na-/qv)dqy at fixed concentrations 

ria/n (or fixed n^/q^). In terms of the inverse matrix element M~^^ the compressibility 

reads: 

18/3M-1 ,. , 

K.T = — —, (A.17 

(8/3e + 9n)Mr/-2/32 ^ ' 

Combining these results one finds the sound velocity as 

c^ = = Tl^ M-^ + n + T15 A.18 

my Cy kt 4 p-^ m,y 9 m,y 6 p m.y 

Alternatively, the right-hand side may be written in terms of the enthalpy per volume 

/i„ = |e + \n/(3. 
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Appendix B. The eigenvalue problem in kinetic theory 

In the main text it has been stated that the frequencies 2™(fc) and their corresponding 
modes of (A;) follow from the eigenvalue problem 1)6. 9|l in arbitrary order of fc. However, 
for our purposes it is more convenient to introduce adjoint modes a™{k). These 
adjoints are particular independent linear combinations of the basis functions a^{k). 
The adjoints follow from the orthonormality condition 

^{[aT{k)raf{k))^S., (B.l) 

With the use of these adjoints, the eigenvalue problem 1().9|) up to second order in the 
wave number may be written as 
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In this expression we have made explicit, through the notation between round brackets 
in the superscripts, that the thermal averages (. . .) have to be evaluated up to second 
order in the wave number. Contrary to what one might expect from the above 
expression, the modes and their adjoints up to first order in k are sufficient to obtain 
the frequencies up to second order. This will be demonstrated below. 

Let us expand the second-order modes a™ (fc) in terms of their first-order 
counterparts a™ (fc)- Formally this expansion can be written as 

a^'^'\k)^Y.^,,af^'^{k) (B.3) 

where the coefficients C^ have the general form 

a, = 6,, + e elf (B.4) 

In the same way we can expand the adjoints: 

C<^)(fc)=5:A,af^^(fc) (B.5) 

j 
with 

A, = % + k^ Ofj (B.6) 

The modes a" (fe) and their adjoints are orthogonal up to first order in the wave 
number, so that one has: 

l{[ar^'\k)raf'\k)r^=&.,+eF^^ (B.7) 

From the orthogonality of the modes a™ (fc) and their adjoints in second order it 

[21 

then follows that the coefhcients Flj obey the relation 

Ff^-{I^T-Cf (B.8) 
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Substitution of the expansions HB.3|I and (|_B.5|I into the eigenvalue equation (|ij.2|l 
yields up to second order in k: 
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Upon approximating the frequencies z™ by z™ in the denominator of the second 
term at the right-hand side, we may write the expression between curly brackets as 
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which defines G\^(. Inserting ljB.10|) in ljB.9p and using the relations (|B.4p . IIB.6|l and 
l|B.8|l . we obtain for i = j: 
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or, with ljB.7p and again (|B.10p 
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(B.12) 



From this expression for the eigenvalue problem it is obvious that the modes and 
their adjoints up to first order in the wave number are sufficient to determine the 
second-order contributions to the frequencies. 
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